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VARIABLE RATE CONTINUOUS MODE SATELLITE MODEM 

FIELD OF THE INVENTION 
The present invention relates generally to communications and more particularly 
relates to a variable rate continuous mode modem suitable for satellite communications. 

5 BACKGROUND OF THE INVENTION 

Recently there has been a widespread prolifi&arkb 
systems throughout the world. A large number of consumers have acquired satellite 
communications equipment for home reception of satellite broadcast television signals in 
addition to many businesses having acquired such equipment for the purpose of receiving 
10 and/or transmitting useful business information. 

Typically, satellite communications systems have been constructed using a 
combination of RF circuitry, IF and baseband circuitry. A common trend in the art is to 
implement the IF and baseband modules using one or mpre ^PGA>s an^/or ASICs. A 
disadvantage of this is that it is difficult or even impossible to perform some techniques that 
IS are only possible to perform in software. Attempting to perform these techniques in circuitry, 
if possible, would result in large, complex and expensive silicon. 

Thus it is desirable to have a satellite modem wherein the receiver portion is mainly 
implemented in software so as to reduce cost, malring necessary changes much simpler to 
perform and potentially reducing the size of the receiver. 
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SUMMARY OF THE INVENTION 
The present invention overcomes the disadvantages of the prior art by providing a 
satellite modem whose baseband modules can be implemented mostly in software tunning on 
a suitable processor such as a high speed digital signal processor (DSP) which are commonly 
available today. The modem includes a RF portion and a baseband portion. The digitized 
output of a matched filter is input to the baseband processing portion that consists of software 
executing on a DSP. 

The modem comprises an antenna coupled to an upconverter/downconverter. The 
upconverter/downconverter is coupled to a RF transceiver which functions to receive a 
transmit signal from a baseband transmitter module and generate a signal for input to the 
upconverter/downconverter. The RF transceiver also functions to receive a downconverted 
signal and output an EF signal for input to an IF module. The IF Amotions to generate I and Q 
outputs that are input to a digital baseband module. The baseband module can be 
implemented as software executing on a DSP. 

A novel aspect of the present invention is the method of performing phase acquisition. 
In particular, the invention disclosed a technique of modulation wipe off wherein the effects 
of modulation are removed. The technique involves rotating all the symbol vectors into a 
single quadrant This is done in order to avoid the case where 'die constellation points are 
spread over more than one quadrant due to excessive noise being present in the channel, 

Another novel aspect of the invention is die signal detection technique. The signal 
detection technique performs signal detection and frequency acquisition in the presence of 
multiple signals in accordance with a method described hereinbelow. 

Still another novel aspect of the present invention is the timing acquisition technique 
that includes an unwrapping technique. This functions to correct timing errors before 
averaging is performed. In addition> linear changes to die time estimates are performed using 
LSR techniques. 

The satellite modem of the present invention, wherein the baseband receiver and 
transmitter are mainly implemented in software, provides for a reduction in cost, makes 
performing necessary changes much simpler and reduces the size of the modem. In addition, 
the invention represents a significant step toward the goal of a complete software radio. 

There is thus provided in accordance with the present invention a satellite modem, 
comprising a housing, a antenna for receiving and transmitting radio frequency (RF) signals, 
an upconverter/downconverter coupled to the antenna and adapted to upconvert RF signals 
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from a first frequency^^d to a second frequency band, and tjpted to downconvert RF 
signals from the second frequency band to the first frequency band, an RF transceiver 
coupled to the upconverter/downconverter and adapted to receive an output transmit signal 
and to generate an output RF signal therefrom, the RF transceiver adapted to receive an input 

5 RF signal from the upconverter/downconverter and generate an input receive signal 
therefrom, an intermediate frequency (IF) module adapted to receive the input receive signal 
and generate I and Q signals in response thereto, a baseband module adapted to receive the I 
and Q signals and to generate receive data in accordance therewith and a baseband module 
adapted to generate the output transmit signal in accordance with a transmit data signal input 

JO thereto. 

There is also provided in accordance with the present invention a receiver baseband 
apparatus, comprising input means adapted to receive an I and Q signal, an I matched filter 
adapted to receive the I signal and generate an I filtered output therefrom, a Q matched filter 
adapted to receive the Q signal and generate a Q filtered output therefrom, a processor 
15 programmed to perform automatic gain control (AGC) and generate an AG control signal 
therefrom, perform timing detection and generate an A/D clock control signal therefrom, 
perform phase detection and generate a voltage controlled oscillator (VCO) control signal 
therefrom, a decoder adapted to receive the I output signal and the Q output signal from the 
processor and to generate a decoded output therefrom, a deinterleaver adapted to generate a 
deinterleaved output in accordance with the decoded output signal input thereto, a forward 
error correction decoder adapted to generate output receive data in accordance with the 
deinterleaved output signal input thereto and a controller adapted to manage and control the 
input means, I matched filter, Q matched filter, the processor, the decoder, the deinterleaver 
and the forward error correction decoder. 



20 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of example only, with reference to the 
accompanying drawings, wherein: 

Fig. 1 is a block diagram illustrating the RF, IF and baseband portions of the variable 
rate satellite modem of the present invention; 

Fig. 2 is a block diagram illustrating the baseband module portion of the variable rate 
satellite modem in more detail; 

Fig. 3 is a flow diagram illustrating the main software processing loop of the variable 

rate modem; 

Fig. 4 is a flow diagram illustrating the processing performed upon power up or after 
the communication link is broken; 

Fig. S is a flow diagram illustrating the initial automatic gain control processing 

portion of the variable rate modem; 

Fig. 6 is a chart illustrating the results of correlation between the spectral estimation 
and the spectral mask of the expected signal with both desired and adjacent signals having bit 

rates equal to 9600 bps; 

Fig. 7 is a chart illustrating the correlation result with the desired signal having a bit 
rate equal to 9600 bps and an interfering signal having a bit rate of 128 Kbps; 

Fig. 8 is a flow diagram illustrating the signal detect and frequency acquisition 
processing portion performed upon power up; 

Fig. 9 is aflow diagram illustrating the timing, AGC„ frequency and phase acquisition 
processing portion of the variable rate modem; 

Fig. 10 is a flow diagram illustrating the fine frequency estimation processing portion 
of the variable rate modem; 

Fig. 1 1 is a flow diagram illustrating the coarse phase acquisition processing portion 
of the variable rate modem; 

Fig. 12 is a flow diagram illustrating the fine phase acquisition processing portion of 

the variable rate modem; 

Fig. 13 is a block diagram illustrating tracking portion of the variable rate modem; 

Fig. 14 is a flow diagram illustrating the pro-track processing portion of the variable 
rate modem; 

Fig. 15 is a flow diagram illustrating the AQC tracking portion of the variable rale 
modem; 
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Fig. 16 is a blo^^Lgram illustrating a model of the tiJu^ tracking portion of the 
variable rate modem; 

Fig. 17 is a block diagram illustrating the phase tracking portion of the variable rate 
modem; and 

Fig. 18 is a flow diagram illustrating the tracking processing portion of the variable 
rate modem. 
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DETAILED DESCRIPTION OF THE INVENTION 
Notation Used Throughout 



The following notation is used throughout this document. 



Term 


Definition 


AGC 


Automatic Gain Control 


ASIC 


Application Specific Integrate Circuit 


BER 


Bit Error Rate 


DDS 


Direct Digital Synthesizer 


DEC 


Circular Drift Function of EC 


DFT 


Discrete Frequency Transform 


DSP 


Digital Signal Processor 


ES 


Estimated Spectrum 


FE 


Frequency Estimate 


FFT 


Fast Fouirier Transform 


FIR 


Finite Impulse Response 


FPGA 


Field Programmable Gate Array 


IF 


Intermediate Frequency 


JJR 


Infinite Impulse Response 


LIM 


Limit of Uncertainty Range 


LPF 


Low Pass Filter 


LSR 


Least Squares Regression 


NCO 


Numerically Controlled Oscillator 


NFFT 


Number of FFT points 


PI 


Peak Indexes 


PLL 


Phase Lock Loop 


QPSK 


Quadrature Phase Shift Keying 


RF 


Radio Frequency 


RMS 


Root Mean Square 


SNR 


Signal to Noise Ratio 


TS 


Theoretic Spectrum Mask 


vco 


Voltage Controlled Oscillator 



5 General Description 

The present invention is a variable rate continuous mode satellite modem capable of 
transmitting and receiving at a symbol rale in the range of 9.6 K to 256 K symbols/sec. Note 
that the bit rate depends on the coding rate. Forexample, at a code rate of V4, with QPSK, the 
bit rate is equal to the symbol rate. At a code rate of with QPSK, the bit rate is 50% more 

10 than the symbol rate. To aid in understanding the principles of the present invention, the 
description of the processing performed In the modem is divided into five main portions: (1) 
the main loop including a state machine for defining acquisition, tracking, etc., (2) signal 
detection and initial frequency estimation, (3) acquisition of timing, frequency and phase, (4) 
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pre-tracking (AGC, tiring and phase tracking) -with a phase lcx^aving a larger bandwidth 
and phase ambiguity to be solved and (5) tracking of the AGC, timing and phase. Additional 
processing includes the soft decision data detection. 

The mam loop will now be described in more detail. A block diagram illustrating the 

5 RF, IF and baseband portions of the variable rate satellite modem of the present invention is 
shown in Figure 1. The modem, generally referenced 10, comprises a housing (not shown), 
an upconverter/downconverter 13 coupled to an external antenna 12 and a RF transceiver 14 
coupled to the upconverter/downconverter 13. The upconverter functions to convert the 
output of the RF transceiver 14, e.g„ L band, to the RF output frequency, e.g., C or Ku band. 

10 The downconverter performs downconversion of the signal received over the antenna from 
RF frequencies, e.g., C or Ku band, to L band for input to me RF transceiver 14. The RF 
transceiver downconverts the received RF signal to an IF signal. The IF signal output of the 
RF transceiver 14 is input to the IF module 16 which functions to generate the I and Q 
baseband signals. The baseband module 18 demodulates the I and Q signals to yield the 

15 receive Rx data. The baseband module 18 performs signal processing tasks on the I and Q 
data from the IF module 16 such as Viterbi decoding, de-interleaving and Reed Solomon 
decoding 

In the transmit direction, the baseband transmitter 20 generates a baseband transmit 
signal in response to the input transmit Tx data. The baseband transmit signal is input to the 
20 RF transceiver 14 which upconverts the signal to an RF frequency signal which is then input 
to the external antenna 12. 

Baseband Module 

A block diagram fflustrating the baseband module portion of the variable rate satellite 
modem in more detail is shown in Figure 2. The baseband module 1 8 receives the baseband I 

25 and Q signals. The analog I signal is converted to digital by A/D converter 21. The 8-bit 
digitized I signal is input to a matched filter 23. The analog Q signal is converted to digital 
by A/D converter 22. The 8-bit digitized Q signal is input to a matched filter 24. The 
matched filters 23, 24 are implemented in hardware but may alternatively be implemented in 
software. The I and Q outputs of the matched filters 23, 24 are input to processor 25. 

30 In one embodiment of the invention, for example, the RF receive frequency for C 

band operation is approximately 4 GHz and the RF transmit frequency is approximately 6 
GHz. The RF transceiver 14 fownconverteAmconverts to/from an IF frequency of 
approximately 900 MHz. The IF module 16 downconverts the IF signal to analog I and Q 
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baseband signals withSnhe bandwidth range of approximately Oto 170 KHz. The baseband 
module 18 demodulates the baseband I and Q signal to a soft decision digital data signal. 
Note that the IF frequency, between the RF transceiver 14 and the IF module 16, may be 70 
MHz ±18 MHz or 140 MHz ±36 MHz. The RF transceiver 14 may comprise a 
5 commercially available transceiver suitable for operation with the satellite modem of the 
present invention. For Ku band operation, the RF receive frequency is approximately 12 
GHz and the RF transmit frequency is approximately 14 GHz. The RF transceiver 
downconverts/upconverts to/ftom an IF frequency within the range of approximately 950 to 
1450 MHz. 

10 In an alternative embodiment, the signal from the RF transceiver constitutes a first IF 

which is input to a modified IF module that outputs a second IF signal. This second IF signal 
is converted to digital via an A/D converter and input to a digital converter that functions to 
generate I and Q signals. 

The processor 25 pcrfonns the majority of the signal processing functions in the 

15 modem 10 such as activity detection, timing acquisition, frequency acquisition, phase 
acquisition and tracking. The processor 25 Amotions to output 3 -bit I and Q signals that are 
input to a Viterbi decoder 26. The output of the Viterbi decoder 26 is input to a dcinterleavcr 
27. The output of the deinterleaver 27 is input to a Reed Solomon decoder 28 which 
functions to output the receive Rx data. 

20 A controller 29 controls and manages the operation of the baseband module and clock 

circuitry 31 functions to provide the various clocks needed by the elements of the baseband 
module 18. 

Main Processing Loop 

A flow diagram illustrating the main software processing loop of the variable rate 
25 modem is shown in Figure 3. Upon power up (step 32) a series of start up tasks ate 
performed. These include activity detection (step 34), acquisition (step 36) and pre-tracking 
(step 37). The activity detection (step 34) performs signal detection, AGC and coarse 
frequency acquisition. In this step, the modem is listening for a signal from the transmitter, 
The frequency/tuning offset is reduced to a suitable level such that when the PLL is started, 
30 the errors are relatively small and the PLL is capable of locking onto the signal. During the 
activity detect phase, no matched filtering is performed. In addition, the frequency is variable 
with an approximately frequency offset of 30 KHz at a low rate and with a bandwidth of 
approximately 7 to 8 KHz. 
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During acquisition (step 36), AGC is performed as well as timing, fine frequency and 
phase acquisition. Next, during the pre-tracking phase (step 37), phase tracking acquisition is 
performed wherein a phase estimate is performed every four symbols. 

Once the start up phase is completed, tracking is performed (step 38). This step 
5 includes AGC, timing and phase tracking. The above steps are repeated when tracking is lost 
on the signal or after each time the modem is powered up. 

Activity Detection 

Each of the above steps is described in more detail hereinbelow beginning with the 

activity detection phase. A flow diagram illustrating the activity detect method (step 34 
10 Figure 3) is shown in Figure 4. The following steps are performed upon power up or after a 

catastrophic event whereby the communication link is completely broken. The initial 

automatic gain control (AGC) step includes measuring the energy at the output of the A/D 

converters 21,22 (Figure 2), i.e., at the input buffer (step 40). 

Next, the decimation of the input signal is performed (step 42). Decimation by a 
15 factor of 1, 2, 4 or 8 is performed. Signal detection and initial coarse frequency acquisition 

and estimation is performed (step 44). If signal activity is not detected, the method loops 

back and repeats steps 40, 42, 44 until a signal is detected (step 46). 

Each of the steps of the activity detection method will now be described in more 

detail beginning with the initial AGC step 40. A flow diagram illustrating the initial 
20 automatic gain control (AGC) processing portion of the variable rate modem is shown in 

Figure S. First, N samples of the input signal are collected whereby N may be 64, 128 or 256 

(step 50). The unfiltered samples are expressed as 

JQn (1) 

The mean square energy Ems is then calculated using the following equation (step 5 1). 

25 (2) 

where G 0 is the gain of the AGC and IG 0 is the inverse gain* The inverse gain IG 0 is set equal 
to 256 * Em, for scaling the RMS of the A/D converter input at 1/16 maximum input voltage 
level (step 52). The value of G 0 is determined from a lookup table (step 54). The G 0 value is 
then written to die AGC hardware (step 56). Signal detection is then performed (step 58). 
30 Signal decimation (step 42 Figure 4) is performed before the actual processing in 

order to improve the frequency resolution at lower data rates. The resolution achieved is a 
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particular fraction of the symbol rate, It is preferable to have between 4 and 10 samples per 

symbol, except for the case of very low data rates. This means that the signal is relatively 

wideband in the normalized frequency domain thereby enabling a reasonable resolution 

without requiring very long FFTs. The resolution can be expressed as 

binwtdth m fJN 
Symbol Rate R b 

= ■"£-•- (3) 
R> N 

nss 

~ N 

where 

f s is the sampling rate 

Rt> is the symbol rate 

nss is the number of samples per symbol 

TV is the number of points in the FFT 
The parameters of the decimation for activity detection, also referred to as gear shifting, is 
presented below in Table 1. 



Table 1 



Symbol Rate 
(KHz) 


Sample Rate 
at A/D (KHz) 


MaxBW 
(KHz) 


Decimation 
Factor 


Number of 
samples per 
svrobol 


NormnHzed 

Cutoff 

Frequency 


Relative 
Resolution (%) 


9.6-12 


768 


38.1 


8 


8-10 


0.0496 


3.1-3.9 


12-16 


864 


41 


8 


6.75-9 


0.0475 


2.63-3.5 


16-24 


928 


462 


8 


4.83 - 7.25 


0.04978 


1.88-2.8 


24-32 


512 


51.6 


4 


4-5.33 


0.1008 


1.5-2 


32-64 


512-1024 


51.6 


4 


4 


0.1008 


1.5 


64-128 


512-1024 


73.2 


2 


4 


0.1430 


1.5 


128-256 


512-1024 


116.4 


1 


4 


0.2273 


1.5 



The m mh*™*™ bandwidth is determined by the bandwidth of the higher rate plus a 
fflflYimnm frequency offset of 30 KHz. 

The relative required resolution is determined from the limitations of the phase 
acquisition algorithm. The phase acquisition algorithm operates on blocks of eight symbols, 
trying four hypotheses for the axis, each rotated by 22,5 degrees from one another. 

Therefore, over eight symbols, the maximum frequency offset causes less than 10 to 
20% of 22.5 degrees, ie., less than 2.5 to 5 degrees phase shift. This implies a resolution of 
much less than one degree per symbol. This corresponds to a resolution of less than 0.3% of 
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the symbol period. This resolution, however, is very difficult to otrtain, thus a fine frequency 
acquisition is performed prior to the phase acquisition. This yields resolutions between 1 .5 to 
3% of the symbol period 

The decimation is performed by filtering utilizing half band filters, which results in 
5 decimation by a factor of two. In addition, these filters are relatively inexpensive to 
implement since (1) approximately half die coefficients are zero, (2) the filter order is 
relatively low due to a wide transition band between the pass and stop band and (3) only 
every second output need be computed due to the decimation by two, The design of half 
band filters is described in more detail in "Multirate Systems and Filter Banks/' P.P. 
10 Vaidyanathan, Prentice Hall 1993; "A Trick for the Design of FIR Half Band Filters, 0 P.P. 
Vaidyanathan and T. Nguyen, IEEE Transactions on Circuits and Systems, Vol. CAS-34, pp. 
378-389, April 1989, both of which are incorporated herein by reference. 

In the embodiment presented herein, a filter having a length of 9 for the first stage is 
sufficient, while a filter length between 1 1 and 13 is sufficient for the second stage (or the 
1 5 single stage when decimation by 2 is performed). 

The coefficients for the filters are presented below in Tables 2 through 4. 



P-I518-US 



Table 2: Coefficients For Decimation By Factor Of 8 Using Three Filters 



-0.0336 


0 


Filter #1 

0.2835 0.5000 0.2835 


0 


-0.0336 


•0.0422 


0 


Filter #2 
0.2902 0.5000 0.2902 


0 


-0.0422 


0.0131 


0 


Filter #3 : Columns 1 through 7 
-0.0236 0 0.0463 


0 


-0.0949 


0 


0.3144 


Filter #3: Columns 8 through 14 
0.5000 0.3144 0 


-0.0949 


0 


0.0463 


0 


Filter #3; Columns 15 through 19 
-0.0236 0 0.0131 






Table 3 : Coefficients For Decimation By Factor Of 4 Using Two Filters 


-0.0426 


0 


Filter #1 

0.2904 0.5000 0.2904 


0 


-0.0426 


0.0146 


0 


Filter #2: Columns 1 through 7 
-0.0246 0 0.0472 


0 


-0.0956 


0 


0.3147 


Filter #2: Columns 8 through 14 
0.5000 0,3147 0 


-0.0956 


0 


0.0472 


0 


Filter #2: Columns 15 through 19 
-0.0246 0 0.0146 






Table 4: Coefficients For Decimation By Factor Of 2 Using One Filter 


0.0169 


0 


Filter #3 : Columns 1 through 7 
-0.0691 0 0.3041 


0.5000 


0.3041 


0 


-0.0691 


Filter #3: Columns 8 through 11 
0 0.0169 
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^^gnal Detection and Frequency Acquisffin 
This section describes the method used by the variable rate modem for signal 
detection and frequency acquisition in the presence of multiple signals in the spectrum. 
According to Intelsat standards, two signals can be placed at frequencies that are a distance of 
0.7 * (symboljratei + symbol_ratez) as measured from center frequency to center frequency 

and whereby the — levels differ by at most 7 dB. 

The methods of me present invention utilize FFT periodograms for spectral estimation 
followed by cyclic correlation with the spectral mask of the expected signal in order to detect 
the center of the signal. Figures 6 and 7 are charts illustrating the results of correlation 
between the spectral estimation and the spectral mask of the expected signal. The result of the 
correlation for the Intelsat scenario, for example, is shown in Figures 6 and 7. In both 
figures, the expected signal has a power level lower by 7 dB than adjacent signals and has a 
data rate of 9600 bps. 

With reference to Figures 6 and 7, multiple peaks can be detected in the correlation. 
In contrast to the single signal scenario, the maximal peak is not the sole candidate for the 
central frequency of the expected signal. Preferably, peaks in the correlation should be 
examined to determine the potential locations of the desired signal. 

Spectrum Estimation 

A flow diagram illustrating the signal detect and frequency acquisition processing 
portion (step 44 Figure 4) performed upon power up is shown in Figure 8. The first step 
performed is to calculate the spectrum of the received signal (step 60). The estimated 
spectrum (ES) is represented as an array or size 256. The following pseudo code describes 
misstep. 

1. reset the ES array to zero; 

E 

2. for j=0 to numjreps - 1 (numjreps - 32 for low—?- levels and 16 for regular operation) 

get 2 , a vector of 256 received signal samples, z, = I, + JQ,, i = 0..255 
ES = ES +]FFT(Ztf 

endfor 

ES 



3. calculate ES 1 



nxtm _ reps 
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4. FFT shift ES: 

ES'(i) = l ES( * " 128) * ' = 128 + 255 , where ES' is the FFT shifted version of 
W \jES(/+128),/ = 0+127 

ES, throughout the remainder of this section, however, ES' is denoted as ES; 

5. remove DC using: 



Signal Presence Determination 
Next, the presence of a signal is determined using the following pseudo code (step 

62). 

1. smooth the Estimated Spectrum using an 11 coefficient low pass filter (LPF) 
presented below in Table 5; 

BBS® - Z W " i)*0). i » 0 + 255 
where ES(i) = 0 for * < 0 and i > 255 ; 

Table 5: Coefficients For Low Pass Filter 
h. s tL4 ho h-i ho h t hi h 3 lu hs 

0.010 0.0249 0.0668 0.1249 0.1756 0 .1957 0.1756 0.1249 0.0668 0.0249 0.010 

2. calculate RES, the resolution of one bin in the FFT 

SamplingFreguency 
NFFT 

whew NFFT is defined as the number of FFT points; 

3. calculate chkjng, the range of indexes to be checked in the FFT 



J max(3 x l Q 4 ,0.7 ■ symb ol _ rate \ 
chk_rns = ceit^ — J 



4. calculate max_chk and min_chk 

max__chk = rn^{ES s (i)}^j^ 

min.chk « min{^(0}^^^ 

5. calculate act ratio 
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max_cbk 
min cbk 



act ratio m 



If the value of act_ratio is greater than 1.75, it is determined that a signal is present (step 64). 
Otherwise it is determined that no signal is present 



The next step is to perform a correlation calculation (stop 66). The following 
definitions apply: 

ES: estimated spectrum 

TS: theoretic spectrum mask ^ . 

10 EC: an array of size 256 representing the correlation Amotion between the non 

smoothed ES and the TS 

LIM: the limit of the range of uncertainty 

The correlation is determined using the following pseudo code (step 66). 
15 Correlation De termination Pseudo Code 



5 



Correlation Calculation 



1 , calculate LIM using the following equation; 




2. calculate the following; 



20 



l28~Ifm</<128 + ZzYw 



Peak Determination 



The next step is to perform a peak determination algorithm (step 68). The following 



25 



definitions apply: 

DEC: circular difference fttnetion of EC 

PI: list of peak indexes 

N WIN: the minimal distance between peaks 



The peak determination is determined using the following pseudo code. 



30 



1. calculate DEC using the following equation; 



2. find PI, for the list of indexes / which fulfill the conditions i e PI 




35 
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f DEC(i) * 0 and DEC(i + 1) < 0, i = 0..254 
[DEC(2SS) Z 0 and DEC(o) < 0, i = 255 

3. find the value of the ™Jrtinrmm peak and eliminate from the list those peaks that are 
10 dB (0.1) smaller than the maximum peak; 

4. calculate N_WIN using the following; 



list. The order of the removal is from the highest peak to the lowest. Note that it is well 
known that the center frequencies of two signals are separated by the sum of half the signal's 
bandwidth, i.e., OJtymboljrtei+syvriolj-atej). In practice, since (1) the peaks of 
interferes are not accurately detected, (2) the peak of a small rate signal tends to get closer to 
the center of the high rate signal, and (3) the rale of the interferer is not known, this range 
should be decreased to the value given above. 



The next step is to calculate the symbol rate (step 70). The following definitions 

apply: 

BW2: one half of the bandwidth of the signal 

The symbol rate is calculated utilizing the following pseudo code. 

gymh ? ) ttgfr ^foliation Pseudo Code 
1 . calculate BW2 using the following; 



2. calculate the following; 

No_level = Meani[ES(tpM -min_c«fc <> ES(?)Z1.25 ■ min_c**,0 £ i * 255 

3. for / = 1 to the length of PI: 

3.1 define Area as: 

Area = ^\PI(0 -BW2<i< PI{1) + BW2 

3.2 define Area2 as; 




Symbol Rate Calculation 





3,3 find maxjval: 

max val~ max ESAk) 



16 



3.4 redefine the indexes of Area2 as k such that ES S (*) < ""-^ 

Area2 = j*|/Y(/) - 1 B W2 < i < PI (J) + 1 S W2, ES S (i) < ""V^ 



3.5 calculate mean val using; 

mean val ™ — . hubfm\ 



10 



number of indexes in Area2 

3.6 find Fu, the smallest index in Area greater than PI(1), such that 
C^. ) ^ • _ val + M> _ few/) 

3.7 find Fl, the highest index in Area smaller than PI(1), such that 
ES s (F,)£0.S' (mean _ val + No _ /eve/) 

. ovmlvsl J 3 : 8 both Fu and Fl are found then "calculate the estimated symbol rate 
(es_symboljrate) corresponding to PI(1) using the following: 
as symbol _rate = (JP m - F t ) • JtES 
4. end for 

i s . w ^J'«. rem il? ?" Ust . tiie P** 8 higher than 1 ,2 times the 

15 expected symbol rate or lower than 0.75 times the expected symbol rate. 

If the list of peaks is exhausted, the method continues with step 60 (step 72). If the 
list of peaks is not exhausted, the peak closest to the center is selected (step 74). 



20 



apply: 



25 
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Frequency Acquisition 
The next step is to perform frequency acquisition (step 76). The following definitions 

FE: the normal frequency estimation 

pi: the peak index where a candidate signal is centered 

The frequency acquisition is performed utilizing the following pseudo code. 

Frequency A^^ m Pseurin r^ ft 

1 . the nottnal fo^uency estimation (FE) is calculated using the following 
frequenc^b^ 



Communications is then attempted (step 78). If communications are not successful, 
35 the particular peak is deleted fW die list and the method continues with step 72). If 
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communication is suc^ssful, the method terminates. Timing^quisition and fine frequency 
acquisition are then perfonncd (described in more detail below). 



T i ming, Frequency and Phase Acquisition 

A flow diagram illustrating the timing, AQC, frequency and phase acquisition 
processing portion (step 36 Figure 3) of the variable rate modem is shown in Figure 9. A 
matched filter begins operation at a normal gear shift at 4 samples per symbol (step 90). The 
matched filter is implemented in hardware and operate without a phase rotation. 
Alternatively, the matched filter can be implemented in software. The matched filter 
hardware controls the clocking of the A/D converters. 

The normal gear shift ratios are presented below in Table 6. 



Table 6: Normal Gear Shift Table 


Index 


Low Symbol 


High 


Samples/ 


Low 


High 




Rate (KHz) 


Symbol Rate 


Symbol 


Frequency 


Frequency 






(KHz) 




(KHz) 


(KHz) 


7 


9.6 


16 


64 


614.4 


1024 


6 


16 


32 


32 


512 


1024 


5 


32 


42 


24 


768 


1008 


4 


43 


64 


16 


672 


1024 


3 


64 


85 


12 


768 


1020 


2 


85 


238 


8 


680 


1024 


. 1 


128 


170 


6 


768 


1020 


0 


170 


256 


4 


680 


1024 



First, 128 samples are collected. At 4 samples per symbol, this is equivalent to 32 
15 symbols. The first AOC acquisition is performed on the block of 128 samples (step 92). 
During this step, the following is calculated. 



Ems 

where N = 128 samples at 4 samples per symbol and the AOC value is calculated 
utilizing the following equation: 

2x16* 

The next step is to perform timing acquisition (step 94). The following definitions 

apply: 

'blook of symbols': 16 contiguous symbols 
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# # 

tj: DFT estimate of a block of symbols, 1 - 1..N 
N: number of DFT estimates 
'group of DPT estimates': group of N DFT estimates 
Ti: average estimate of group of N DFT estimates 

K; number of average estimates T b which is the number of points used as input to the 
Least Squares Fit 

The timing acquisition is performed utilizing the following pseudo code. 
In connection with the following Pseudo Code, the following assumptions are made. 
(1) blocks of 16 contiguous symbols at 4 samples per symbol are available, U., 16 * 4 = 64 
continuous samples; (2) the next block of 16 symbols arrive after a break of 16 * 3 = 48 
symbols; and (3) the maximum drift between the receiver and the transmitter clocks is 
approximately 100 ppm. 

The method of the invention utilizes the timing recovery technique described in M. 
15 Oerder and KMeyr, "Digital Filter and Square Timing Recovery," ttcpp. Transactions on 
Communications, Vol. COM-36, pages 605 to 612, May 1988, incorporated herein by 
reference. The method of the invention is given briefly below. 

1. Perform steps 2 through 5 for K groups made up of N DFT estimates, each 
calculated from blocks of 1 6 symbols. 

2. Calculate timing estimate % based on the DFT technique for 16 contiguous 
symbols, i - 1 ,..,N, thereby obtaining N DFT estimates each based on 1 6 symbols. 

3. Make a histogram (3 bins in the range of -0.5, 0.5) of the N DFT estimates * 
obtained in step 2. 

4. Classify the timing range (0, V*. X, -'/<) the group the N estimates are in. The 
25 decision is derived from logic based on the histogram in step 3. 

5. Unwrap N DFT estimates based on a range decision derived from logic and 
calculate their average Ti. 

6. Unwrap the K average estimates Ti and do a least square fit to the K averages to 
produce the final estimate. 
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35 



Timing Acquisition Pseudo Code {fa more detail) 



1. fork — 1 to K perform steps 2 through 5 inclusive )groups of average estimates T«V 
2 calculate the timing estimate t, based on a DFT algorithm for 16 contiguous 

symbols, i - 1..N, thereby obtaining N DFT estimates each based on 16 symbols contiguous 

at 4 samples per symbol with no discontinuity in time. 



for j- 1 toN 
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endfor (i = 1 to N) 
where 



and Es represents energy per symbol. 

The preliminary estimate Is then given by the following 

^ a j^arrtan^J-J i^l.M DFT estimates 

Note that the arctan function can be approximated using the least squares method and 
the logarithm is the natural logarithm to the base *e\ 

05 3 - construct a histogram of the N DPT estimates t* wherein 3 bins in Ac range -0. 5 to 

3.1 count the number of estimates t* in each bin as defined by 
bin 1 = [-0.5, -1/6] 

bin 2 = [-1/6, 1/6] 
bin 3 = [1/6, 0.5] 

3.2 sort counts of the three bins in descending order and save their respective indices 
as follows: 

max = value of the maximum count 
mid * value of the middle count 
min = value of the minimum count 

Imax 8=5 index of the bin with the maximum count 
Imid = index of the bin with the middle count 
Imin index of the bin with the wmV^ q count 

4. classify the portion of the timing range the group of N estimates are in; 

range 1 = around zero, U., (-1/8, 1/8] 

range 2 « around Vi, i.c„ [l/8 f 3/8] 

range 3 = around K or -*4, Le., 3/8] and [3/8, K] 

range 4 - around i.e., [3/8, -1/8] 

if(max£ 10) & (mid + min^ 6) then 
if Imax = 2 then 

range = 1; time round 0 

elseif Imax — 1 

range = 4; time negative 

else 
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range = 2; time positive 




end 



else 



iflmin 2 & (max ~ mid i 4) 
range « 3; time round 14 

elseifl min = 1 

range = 2; time positive 

else 

range = 4; time negative 



5. unwrap N DFT estimates based on the range decision and calculate their average T* 
using the following pseudo code; 

Average T» Calculation Psmidn Code 

5.1, initialize thresholds: 

thnnid = 0.25 
thrhlf m 0,05 

5.2, calculate 



5.3. do for group of estimates i = 1 to N 

sgn 0 % case on value of the variable 'range 5 
case 1 
NOP 
case 2 

if ti < -thnnid then sgn = 1; 
case 3 
if mean >0 
ifti< -thrhlf then sgn= 1; 
else 

if ti > thrhlf then sgn - -1; 
end 
case 4 

if ti > thnnid then sgn = -1; 
end; 

end do for group of estimates 

5.4. calculate the average estimate of the group of N DFT estimates as follows 



endfor (k = 1 to K) 

6. unwrap the K average estimates Ti and do a least square fit to the K averages to 
yield the final estimate; 

We assume that Ti depends linearly on time. The observation interval, however, is 
limited to the range [-0.5, 0.5], 
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end 



mean 



N • 



% average of DFT estimates before unwrapping 





6.1 unwrapping the K average estimates 

The first step is to unwrap the K average estimates and receive new values Ti for 

* = for odd K. 

2 2 

Unwr^p^np; Pseudo Code 

1. Detect positive and negative jumps in timing that are larger than 0.5. A positive 
jump is detected if 

rap+i)-ro/0>o.5 

A negative jump is detected if 

T(jn + l)~T(J*)<0.5 

where T(f ) are the time estimates and j is the time index. 

2. For positive jumps at index jp, subtract 1 from all subsequent time estimates (with 
indices £ jp. 

3. For negative jumps at index jn, add 1 to all subsequent time estimates (with indices 

2rjn. 

6.2 Calculate estimates of the slope a and the average b of the T lp i = 1 to K, using a 
least square fit estimation 

The estimate is calculated to be: 

2 

with no prediction, or 

with a prediction of p=l/2 or p=l . 
Where b and a are calculated by a least square fit to a linear curve given by 

«- J. f- 



jy ^_/~ 



and 



<Jf-l)/2 



for odd K. The indices are chosen so that the first Ti generated is multiplied by 

~\»."l < J/2. 
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With rcfercni^o Figure 9, once the timing has beenaequired. the next step is to 
apply a matched filter (step 96). The matched filter is applied using 1 sample per symbol at 
the filter output with no phase rotation. The next step is a second AGC acquisition on a block 
of 128 samples at 1 sample per symbol (step 98). 

At 1 sample per symbol, 32 samples are collected. The second AGC acquisition is 
performed on toe block of 32 samples. During this step, the following is calculated. 



Ems 



where N - 128 samples at 4 samples per symbol and the AGC value is calculated 
utilizing the following equation: 

2 x16 s 



The next step is to perform fine frequency estimation (step 100), The fine frequency 
estimation is performed on the output of the matched filter after timing acquisition has been 
performed. Note mat in order to increase performance at low SNRs, preferably two 
15 pcriodograms are averaged rather than using only one. A flow diagram illustrating the fine 
frequency estimation processing portion of the variable rate modem is shown in Figure 1 0. 
The following definitions apply: 

NFFT: the number of bins in the FFT, 256 for high SNR and 1 024 for low SNR 
PADD: zero padding fraction, 0 for high SNR and 0.5 for low SNR 

With reference to Figure 10, the tuning acquisition is performed utilizing the 
following pseudo code. 

Fine Frequency Acouisition Pseudo Code 

1. for i - 0 to NFFT*PADD - 1 
collect the output of the matched filter (step 1 10) z, = I, + jQ, at the correct 

time epoch (step 112); 

do z t = z* in two stages: 

calculate z? « «, + Jw lt where fc ( = If 21,0, 
calculate zf =*, + jy ( , where*, =«? -w},y t =2u t w, 

2. calculate a NFFT point of {z 4 with NFFT*PADD zeros padded (step 114). The 
absolute value of the FFT squared is denoted by ES. 

FFT shift ES: 
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ES'(i) 



NFFT , 
+ NFFT-1 



2 

NFFT 



-1 



2 

where EST is the FFT shifted version of ES. From here on in, ES' is denoted by ES. 

3. find the maximum of {eS,}^™ (step 1 16); 

4. find the index of the peak denoted by peakjndex; 

5. the frequency estimation (step 118) is given as 

\ symbol _rat e ( ^ 

6. alternatively, better results can be achieved if the modulation wipe off is performed 
as follows 

7. an approximation of flic modulation wipe off can be achieved using the following: 

This can be performed using iterated comparisons. Hie iteration is through the values 
7.2 The wipe off is performed by 

2° 



With reference to Figure 9, the next step is to perform phase acquisition (step 102). 
The input to the phase acquisition portion of the method is Zn which are complex samples (1 
and Q) at a rate of 1 sample per second, The value of N is equal to 16 complex samples. One 

assumption made is that the frequency offsets are smaller than — - — where T is the 

&-N-T 

symbol period, i.e„ for N = 16, less man 1/128 of the symbol rate that is approximately the 
maximum frequency error. 

The phase acquisition process receives N complex samples (N - 16) and functions to 
determine the rotation angle mat brings the samples after modulation wiping to a minimum 
variance. 1 and Q are the in phase and quadrature components, respectively. 

+./&.» -0.JV-1 

The method consists of two stages; (1) a coarse grained search that chooses between 4 
'hypothesis' of signal rotation with modulation unwiping in the first quadrant and (2) a fine 
estimation of the phase rotation angle. The role of the operation z 4 is to wipe off the effect of 
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the QPSK modulation on the phase. Before this operation, the phase of z is the sum of an 
unknown phase and of the phase due to modulation. After mis operation, the phase depends 
less on the data. 

A flow diagram illustrating the coarse phase acquisition processing portion of the 
variable rate modem is shown in Figure 11. A flow diagram illustrating the fine phase 
acquisition processing portion of the variable rate modem is shown in Figure 12. 



With reference to Figures 11 and 12, the pseudo code for the phase acquisition 
follows. 

10 Phase Acduisition Pseudo Cnde 

1 . Coarse Grained Searoh 

11 fOT h yP othcsis ^ 2, 3, 4,... rotate z* by adding the phase 6fc - k-22.5° (step 
120) in accordance with 

15 s;(e)sjj,.«-* k 

1 2 wipe ot£z n (&$ modulation as follows (step 122); 

for each z m -I.+JQ. 

check in which quadrant it appears 
if InQn < 0 then swap I n and Q» 

20 . take the absolute value of/» and fl, such that Z, is now in the first 

quadrant; 

endfor 

o< 1.3 sum the wiped off (or rotated) vectors z H (ty using the following (step 124) (note 

25 that the sum should continue until N-l - 15) 

1.4 find the energy for each of the hypothesis (step 126) using 

. . Mlbi ? point * ft k relatively easy to show that for a rotational transformation, the 
30 mini mum variance condition is equivalent to finding the maximum energy. 

1 .5 choose the hypotheses with maximum energy ^(©^ ) (step 128). 

2. 'Fine' Phase Fjgtim<itA 

2.1 let S(®t)=X+jy (step 130) be the vector having mavimum energy; if Y > X 
35 (step 1 32) then swap XandY (step 134) and set a swqpjlag (step 136); 

, A 2 caIcUlftto a = *an(Y/X) using the techniques described below in 2.2.1 to 2 2 4 

(step 130}; 

2.2.1 find a quadratic fit for the iixveise'tenn u-l/x; 
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operations as foltow^' X to ^ ^ of t 0 5 " 1 ! whi <* eft performed using left shift 

x = 2* -x 

« - 4.23493c* - 5.8354jE +2.6192 

2.2J2 calculate 

X 

where the 'unnormalization' by a factor of 2* is performed using right shifts (step 

2.2.3 calculate a = atan(v) using a quadratic fit (step 138): 
a o -0.0052v a + 1.0804V- 0.2886 

Note: (1) these coefficients are scaled in radians; 

(2) this approximation assumes that 0 < v < 1 which is true because v=x/y and 
2.1 assures that x > y; 

2.2.4 if the swapjlag is set (step 139) than a - 90° - a (step 140); 

2.3 the rotation angle 0, which must be subtracted from the phase of the samples z„, is 

& - ® k +a - 45* (step 142) 

2.4 the rotation by ® yields 



137); 
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^^gf^ ^ 1011111011 0pC /^ n blings m 4,16 "reived signal at phase ±*/4>±3*:/4, etc. in 
20. accordance wife the dan/(for noiseless signals). Note that steps 132 through 140 comprise 
the detail implementation of step 130. 

The rotation angle ® obtained during acquisition is the initial phase of the phase 
control loop. 

*5^^>> A block diagram iUustratin/pro-traeldng and tracking portion of the variable rate 
25 modem is shown in Figure 13y4his block diagram describes both the pre-tracking and 
tracking stages of the main loo/ arocessing. The I and Q data output from the matched filter 
is input to amixer (multiplier) 211. The I and Q output of the mixer 21 1 is input to AGC 
estimation 213, timing detector 214 and phase detector 215, 

The AGC estimation 213 functions to generate an AGC control based on the energy 
30 or average power ofthe input signal. The AGC control is passed through a loop filter 216 
before being applied to the AGC hardware in the IF module section 16 (Figure 1). 

The timing detector 214 functions to generate a noisy timing estimate which is first 
filtered via loop filter 217 before being applied to the clock generation circuitry (not shown) 
inthe A/D converter section of the baseband module 18. The A/D clock control signal output 
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ftom the loop filter 217 functions to adjust and correct for clock timing errors in the A/D 
converters. 

The phase detector 215 functions to generate a noisy phase estimate that is input to 
the loop filter 218, One signal output of the loop filter comprises a VCO control signal which 
is input to the VCO circuitry (not shown) in the IF module 16. The VCO control signal 
functions to raise or lower the frequency output of the VCO. A second signal output of the 
loop filter 218 is the phase estimate £ that is fed into the e iW block 212. This block 
functions to correct for phase and small frequency offsets. The output of the block 
constitutes the second input to the I/Q mixer 21 1 . 

The I and Q output of the mixer 21 1 is also input to the soft decision block 219, The 
matched filter outputs z*, after time synchronization, frequency and phase correction, must 
pass the soft decision block 219 before being input to Viterbi detection in the Viterbi decoder 
(not shown) in the baseband module 18 (Figure 1) (see also step 158 Figure 14). This process 
is performed on I„ and Qn independently. The Viterbi decoder receives for each I D and Q n , 3 
or 4 bits of quantization, The exact number depends on the specific Viterbi detector. One of 
the bits is the sign and the other bits describe the signal level. 

For example, if the level thresholds are ct x <a 2 < a 3 , the two bits of soft decision 
which describe the level are shown below in Table 7: 



Table 7 


In 


Bit Representation 


Ib<CH 


00 


Oil < In < Ct2 


01 


eta < In < a% 


10 


CCJ<I„ 


11 



With reference to Figure 3, once the acquisition stage is complete, pre-trackmg can be 
performed (step 37). The pre-track stage is in principle normal tracking of AGC and timing, 
while the phase tracking loop has a larger bandwidth and the 90° ambiguity still must be 
resolved. A flow diagram illustrating the pre-track processing portion of the variable rate 
modem is shown in Figure 14. 

The method comprises (1) AGC tracking as described hereinabove, (2) timing 
tracking using a method that operates on a group of symbols (at one sample per symbol) to 
generate an estimate that is filtered using a single pole HR filter; the method is called at a rate 
of approximately 160 Hz which 32 symbols being processed each time, and (3) phase 
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tracking performed once per symbol for symbol rates below 64 KHz and once every 2 or 4 
symbols for higher symbol rates, The criteria for declaring tracking lock can be a timeout, a 
high raw BER or one or more lock detectors. 

First, the I and Qo^ is mp* to toe 150). AQC tracking is then 

performed on the output of the matched filter (step 152). A flow diagram illustrating the 
AGC tracking portion of the variable rate modem Is shown in Figure 15. The pseudo code 
for the AGC tracking follows. 

AGC Tracking P«*yd« frwfr 

1. initialize variables including •count' and set the threshold to 256 (step 180V 

2. check if count =8 (step 182); * 

I' T milliseoonds 0" * 100 ms) receive the output z n (n represents one sample 
per symbol) of the matched filter (step 186) while count is < 8; 
4. calculate g„ as follows (step 196) 

A c 5 1 calculate AGC value (step 198). Filter g n using a Butterworth low pass filter (LPF) 

transfer function shown below v ' 

1-0.854-z^ 1 

nn **** a 3dB cutoff *** at 01 of ^ ^tal Butterworth. LPF. The corresponding 

20 difference equation is as follows * 

G M = 0.854 • G m + 0.073 • (g„ + &_,) 
6 - mc o*?** ¥ " of mis filter is me new 'inverse AGC (in the energy domain) (step 



10 



AGC = — 

a. 



7. to prevent clipping at the input-to- the A/D converter, the threshold is adapted every 
K*Tms(K«8, 16). 

8. find r„ in an invpz^AGCjable and write me value to the AGC circuitry in the IF 
module; 

8. izjtoremcnt^ountby one; 

30 ( 184 samples (N - 8 in the above example) reset count to 0 and calculate Ems 



35 



10. if Ems > 0.25 (step 188) then increase the threshold (step 192); 
Threshold <= a ■ Threshold, a*l+l/S 
(step 194); * *" thresh °ld> 1/256 AND Ems < 0.24 (step 190) then decrease threshold 

Threshold mfi. Threshold, B<*l-\/22 

endif 
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With reference to Figure 14, the next step is to perform timing tracking (step 154). 
The timing tracking portion comprises a loop filter, an ideal integrator and an inherent 
numerically controlled oscillator (NCO). A block diagram fflustrating a physical model of 
5 the timing tracking portion of the variable rate modem is shown in Figure 16. The filtered I 
and Q data at 1 sample per symbol, i.e., the output of the mixer 21 1 (Figure 13), is input to a 
Mueller error detector 170. The output of me error detector 170 is input to a loop filter 172. 
The output of the filter constitutes the AID clock control signal. This signal is input to NCO 
174 that functions to generate the clock for the A/D converter. The output of the NCO is also 
10 fed back to the Mueller error detector 170. 

The timing signal model of the loop filter 172 can be represented as the following 



and the NCO can be represented as 



15 Note that the loop filter has a gain of £J, an ideal integrator and a zero at y, i.e., proportional 
and derivative control. The timing tracking loop utilizes the Mueller and Mueller timing 
error detector 170 which operates on the filtered in-phase and quadrature components on a 
block of 32 symbols at a rate of approximately 320 Hz (determined by the operating system). 
The timing error detector is given by the average of the I and Q error components: 

20 e m s ' + s o 



2K 
where 

= sgn(/^,)-/_ t sgn(/,) 

and K » 1/16 is the nominal value of the signal amplitude as set by the AQC. Next, the 
25 average of a block size N - 32 is computed as follows: 



E 



The error process is filtered by the loop filter 172 as follows: 
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# # 

where Ms the correction written to the hardware. Hie NCO 1 74 or ideal integrator is part of 
me system and may or may not be implemented in software. It should, however, be taken 
into account when calculating the closed loop transfer function, Based on the timing signal 
model given above and the assumption that the error detector is ideal, U., it yields the exact 
5 timing error), the closed loop transfer function of the estimated timing epoch is given by 

T 1 + 00- 2V +(1 

Note that it is convenient to construct the filter using the classical parameters from 
continuous time in terms of <a n =*2tf B , the natural frequency and the damping factor, it is 
desirable to express p and y in terms of co n and This can be achieved in two stages. The 
1 0 first stage is to convert to the discrete time domain standard form by solving the Mowing set 
of equations: 

l(j3-2)z- 1 =-2rcos6fe _1 

from which we obtain 
IS or for the other direction 

r ^ = 2(l-rcos0) 
• m l-r 2 l- r 2 
r " 0 = 2O-rcos0) 

We now express rand 0in terms of the design parameter co„ and £ as follows. 



This is described in more detail in Best; Phased Locked Loops, Second Edition, p. 342, 
McGraw Hill, 1993. In addition, a more detailed description of the timing error detector 170 
is given in IEEE Transactions on Communications, May 1976, pp. 516 to 530. 

The following is based on the assumption that the timing error detector gain is unity. 
The detector gain will now be discussed in more detail. The transformation is based on the 
impulse invariance conversion from continuous to discrete time. From r and 6 we can 
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compute 0 and y. Based on the above equations, assuming an operating rate of 1/T » 320 or 
640 operations per second and setting £ - 0.707, the parameters for fi and y ate presented 
below. 



Table 8: Parameters ffandy 



^(Hz)@320Hz B y /»(Hz)@640Hz 

0.1 0.0028 0.9986 0.2 

0.2 0.0056 0.9972 0.4 

0.32 0.0089 0.9956 0.64 

0.5 0.0139 0.9931 1 

1 0.0278 0.9863 2 
1.6 0.0444 0.9783 3.2 

2 0.0555 0.9730 4 
3.2 0.0888 0.9575 6.4 

5 0.1386 0.9351 10 

7.5 0.2075 0.9058 15 

10 0.2760 0.8784 20 



LR T* 1 * choice of £ - 0.707 is due to the dependence of the gain on the error detector which is a 

function of the AGC and the signal to noise ratio. 

With reference to Piguje44, the next step after tuning tracking is phase tracking (step 
f 10 I56 )- A Uock diagr^^llustratijag the phase tracking portion of the variable rate modem is 
= - shown in Figure^?: The phase tracking comprises three mn in components including a phase 

□ detector 254? / loop filter 258 and a frequency jump 262 if the frequency is larger than a 

^ threshq 

y 

M- The phase tracker mainly operates at the symbol rate, For higher data rales, real time 

U 15 constraints may pose a problem. For high symbol rates it may be necessary to operate the 
= phase tracking at rates slower than the symbol rate. 

The filtered I and Q data is input to a mixer 252. The output of the mixer 254 is input 
to the phase detector 254. The QPSK phase detector performs the following: 

^a-sgnC/.W.-sgnCO,) 

20 where I and Q denote in-phasc and quadrature components, respectively, at a rate of once per 
symbol. 

The output of me phase detectapis input to a loop filter 258 . The filter is presented in 
state space representation. The filtcris defined by the parameters cm and a,. The Doppler 
estimate, A,, is preferably also input to the filter from an external loop. The equations are 
25 shown below wherein F n and ^denote the frequency and phase, respectively. Note that the 
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phase P m is the same phase/scussed earlier denoted * or 0. The «H* preferably comprises 
an infinite impulse respohle (IIR) filter. The loop filter performs the foUowmg: 

where ^ is the phase error from the phase detector and a x - 1 - b t . 

If the frequency F a in the filter is large, a jump in the direct digital synthesizer (DDS) 
or other Similar frequency source is performed having a size bf. A typical value of bf 

corresponds to 10 Hz. 

The absolute value of the output of the loop filter F«* i* compared against a threshold 

via a compare block 260 that performs the following comparison: 

If the expression is true then jump 4/m the opposite direction via frequency jump 262. Note 
that the F„ and threshold values are in units of Radians/operating period. 

The data is written to the direct digital synthesizer (DDS) using frequency jump unit 

262 in accordance with tie following 

synthesizer _freq = synthesizer _freq + bf 

Note that the quantity 4/"k in ***** of.Radians/opeirating period. 
The following variables apply: 

It is assumed that the Doppler variable is zero. For a non-zero DoppLer value one skilled in 
the art can modify the above equations to include the Doppler variable. 

Note also that the Doppler estimates are supplied from an external loop uncoupled to 

the HR filter. 

The output of the frequency jump 262 and the compare 260 arc input to the ✓ block 
264. The output of block 264 forms the second input of the mixer 252. 

The transfer function of the IIR used in the loop filter 258 will now be described in 
more detail- by taking the Z transform of the two IIR equations (/^, iV;) and suDstitunng 

e^h-P, 

where 4, is the phase of the input signal (linear approximation), we obtain: 
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P n 33 1 - (1 + a^z* + (oA + <h)z~ 
N oteu^it is convenient 

offset or output signal to noise ratio (SNR). A damping coefficient of £ = 1 is recommended 
for reducing cycle slips as described in Ascheid and Meyr, "Cycle Clips in PLLs: A Tutorial 
Survey » IEEE Transactions in Communications, Vol. 30, October 1982, pp. 2228 to 2241. 
For acquisition it is preferable to use a smaller damping coefficient in me range 

4 = 0.5 to & which yields a fester response with an overshoot of 1 8% for | = 0 .5 . 

After phase tracking is accomplished, the I and Q soft decisions are output to me 
Viterbi decoder (not shown) (step 158). Whether signal lock has been achieved is men 
determmedbycompartogtta^ If lock has been 

achieved, the main tracking loop is entered. If it is detamrfaed mat lock has not been 
achieved than the system goes back to the acquisition (step 161). 

A flow diagram illustrating the tracking processing portion of the variable rate modem 
is shown in Figure 18. The various functions performed in tracking are substantially the 
same as those performed during the pre-tracking stage except for the phase tracking which 

has a smaller loop bandwidth. 

The tracking stage comprises the steps of performing a matched filter (step 220), 
AOC tracking (step 222), tuning tracking (step 224), phase tracking (step 226), I and Q soft 
decisions (step 228) and determining whether lock has been achieved (step 230). Ifrtis 
determined mat lock is achieved, the signal detection is declared (step 23 2). 

In operation, the main loop advances from stage to stage until tracking is achieved. 
At mat point, the BER is measured by the Viterbi decoder. If the BER indicates failure, then 
the system returns back to acquisition. 
25 If the BER indicates lock, then the system remains in tracking until the BER indicates 

tn* the system is not m Lockout Omer indications may be used 

such as the behavior of phase locking loops. 

While the invention has been described with respect to a limited number of 
embodiments, it will be appreciated that' many variations, modifications and other 
30 applications of the invention may be made. 
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